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Abstract  .We  study  the  maximum  temperature  rise  induced  by  a  rotating  or  dithering  Gaussian  laser 
beam  on  a  two-layer  body.  First,  we  derive  semi- analytical  solutions  by  solving  the  transient  three- 
dimensional  heat  equation  in  a  semi-infinite  domain  with  insulating  surface.  Second,  we  provide  numerical 
solutions  for  a  two-layer  structure  in  a  finite  domain.  Our  results  show  that  the  maximum  temperature 
rise  can  be  reduced  significantly  either  by  applying  a  thin  layer  of  coating  on  the  protected  object  or  by 
moving  the  beam  or  equivalently  the  object. 
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1.  Introduction 

Laser- induced  heating  plays  an  important  role  in  materials  processing  [11].  The  theoret¬ 
ical  modeling  of  temperature  profiles  induced  by  laser  radiation  in  solids  has  been  exten¬ 
sively  studied.  For  example,  in  [7]  Lax  modeled  the  spatial  distribution  of  the  temperature 
rise  induced  by  a  stationary  Gaussian  laser  beam  in  a  solid  sample  using  a  one-dimensional 


*  Corresponding  author 

E-mail  address:  hzhou@nps.edu(H.Zhou) 

Received  December  20,  2011 


339 


Form  Approved 
OMB  No.  0704-0188 


Report  Documentation  Page 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2012 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

4.  TITLE  AND  SUBTITLE 

Temperature  Rise  in  a  Two-Layer  Structure  Induced  by  a  Rotating  or 
Dithering  Laser  Beam 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Department  of  Applied 

Mathematics, Monterey, CA, 93943 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Abstract. We  study  the  maximum  temperature  rise  induced  by  a  rotating  or  dithering  Gaussian  laser  beam 
on  a  two-layer  body.  First,  we  derive  semi-analytical  solutions  by  solving  the  transient  three-  dimensional 
heat  equation  in  a  semi-in  nite  domain  with  insulating  surface.  Second,  we  provide  numerical  solutions  for 
a  two-layer  structure  in  a  nite  domain.  Our  results  show  that  the  maximum  temperature  rise  can  be 
reduced  signi  cantly  either  by  applying  a  thin  layer  of  coating  on  the  protected  object  or  by  moving  the 
beam  or  equivalently  the  object. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

21 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


340 


HONG  ZHOU* 


integral  based  on  thermal  conduction  in  the  material.  Since  the  closed-form  solution  in  [7] 
is  a  steady  state  solution  to  the  heat  equation,  it  is  limited  to  stationary  or  slow  moving 
beams.  Later  Lax  extended  his  analysis  in  [7]  to  allow  temperature-dependent  thermal 
conductivity  [8] .  Further  theoretical  and  numerical  developments  include  time-dependent 
solution  for  a  scanning  Gaussian  beam  [6],  laser- melted  front  [3],  continuous  wave  (cw) 
laser  annealing  of  heterogeneous  multilayer  structures  [5],  temperature  profiles  induced 
by  a  moving  cw  elliptical  laser  beam  with  the  inclusion  of  the  temperature-dependent 
surface  reflectivity  and  thermal  diffusivity  [9] ,  a  general  analytic  solution  for  the  temper¬ 
ature  rise  produced  by  scanning  Gaussian  laser  beams  [12],  the  Green’s  function  solution 
to  the  heat  equation  for  a  two-layer  structure  with  scanning  circular  energy  beams  [4], 
the  transient  three-dimensional  analytical  solution  of  the  temperature  distribution  in  a 
finite  solid  when  heated  by  a  moving  heat  source  [1] ,  and  the  laser  forming  of  plates  using 
rotating  or  dithering  Gaussian  beams  [13]. 

However,  the  main  purpose  of  this  paper  is  to  find  ways  to  reduce  the  maximum  tem¬ 
perature  rise  induced  by  a  laser  beam,  which  is  relevant  to  military  applications.  In  [15] 
we  have  derived  the  analytical  solution  for  the  temperature  rise  induced  by  a  rotating  or 
dithering  laser  beam  in  a  semi-infinite  domain  by  solving  a  transient  three-dimensional 
non-homogeneous  heat  equation.  In  [14]  we  have  carried  out  numerical  simulations  to 
obtain  the  temperature  rise  induced  by  a  rotating  or  dithering  laser  beam  on  a  finite 
body.  These  studies  have  shown  that  the  maximum  temperature  rise  can  be  reduced  by 
increasing  the  frequency  of  the  rotating  or  dithering  beam.  In  this  paper  we  would  like 
to  extend  our  previous  studies  to  a  two-layer  structure.  In  particular,  we  would  like  to 
investigate  the  effect  of  a  thin  layer  of  coating  on  the  temperature  rise  induced  by  a  laser 
beam. 

We  organize  our  presentation  as  follows.  In  Section  2  we  derive  semi-analytical  solutions 
of  the  transient  temperature  distributions  induced  by  a  rotating  or  dithering  laser  beam 
in  a  semi-infinite  domain.  We  present  our  numerical  results  for  a  finite  domain  in  Section 
3.  Our  main  results  are  summarized  in  Section  4. 

2.  Semi-analytical  solution  in  a  semi-infinite  domain 
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Consider  a  laser  beam  along  the  2  direction,  hitting  a  solid  surface  in  the  x-y  plane  and 
rotating  in  the  x-y  plane.  The  solid  has  a  two-layer  structure  which  consists  of  a  film  of 
thickness  a  in  the  ^-direction  (region  1)  and  a  substrate  of  infinite  thickness  (region  2). 
The  material  of  the  solid  is  assumed  to  be  isotropic  and  the  solid  is  infinite  in  the  x  and  y 
directions.  Figure  1  depicts  the  geometries  of  the  problem  to  be  solved  and  the  coordinate 
system  that  we  have  chosen.  A  semi-infinite  geometry  is  a  reasonably  good  approximation 
if  the  laser  beam  is  small  compared  to  the  object.  Heat  loss  from  the  target  surface  is 
assumed  to  be  negligible  compared  to  conduction  into  the  solid. 


Figure  1.  A  schematic  diagram  shows  a  rotating  laser  beam  on  a  two-layer  object. 


In  general,  the  heat  equation  for  an  isotropic  medium  with  a  temperature-dependent 
thermal  conductivity  K  can  be  written  as 


Q^pz^t) 

PC 


(i) 


f  -ppV-WUW) 
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where  T(x,  y,  z,  t )  is  the  solid’s  temperature  rise  above  ambient  (i.e.  the  difference  between 
the  solid  temperature  and  the  ambient  temperature)  and  it  is  a  function  of  space  (x,  y,  z ) 
and  time  t,  V  is  the  gradient  operator,  p  is  the  density,  C  is  the  heat  capacity,  K  is  the 
thermal  conductivity,  and  Q(x,  y,  z ,  t)  is  the  incident  heat  source  term.  The  temperature- 
dependent  thermal  conductivity  K(T )  can  be  removed  from  the  heat  equation  by  applying 
a  Kirchhoff  transformation. 

The  Kirchhoff  transformation  introduces  a  linearized  temperature  rise  6  defined  as 

,  ,  „„  n/rT, ,  fT  K(T ')  ,  d6  ATT) 

(2)  0(T)  -  9(T„)  +  ^  jr^dT  .or  —  -  , 

where  9(Tq)  and  K(T0 )  are  constants.  Note  that  when  the  thermal  conductivity  is  con¬ 
stant,  then  (2)  implies  that  9(T)  =  T  if  we  choose  9 (To)  =  T0.  Applying  the  chain  rule, 


we  have 


d9  d9  dT  K(T)  dT  dT  K(T0 )  d9 


dt  dT  dt  AT  (T0)  dt 


dt  K(T )  dt 


_  AT(T) 
“  dT  VT  ~  a: (To) 


VT,  or  VT 


A'(T) 


Substituting  the  expressions  (3)  and  (4)  into  the  heat  equation  (1),  we  can  express  the 
heat  equation  (1)  in  terms  of  the  linearized  temperature  rise  9  as 


Q(x,  y,  z,  t)D 
K(T0)  ’ 


where  D  =  K(9)/pC(9)  is  the  thermal  diffusivity  of  the  material,  and  A  is  the  Laplacian 
operator.  Under  the  assumption  of  constant  thermal  diffusivity,  equation  (5)  becomes  a 
partial  differential  equation  with  constant  coefficients. 

For  the  two-layer  structure  depicted  in  Figure  1,  the  laser  beam  hits  the  surface  z  —  0. 
The  linear  heat  equation  for  region  1  is 

96 1  7-)  Afl  Q(x,y,z,t)D1 

<6)  ~m  “ DlAe 1  =  avTo)  ’ 


and  the  linear  heat  equation  for  region  2  is 


-  Do  A9o 


Qip^zD^Dj 

K2(To)  ' 
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where  D i,  D2  are  the  thermal  diffusivities  of  the  materials  in  region  1  and  region  2, 
respectively.  Initially  the  temperature  rise  is  assumed  to  be  zero.  The  boundary  conditions 
for  6i  and  6b  are 


(8) 

(9) 

(10) 
(11) 


dOl 

— —  =  0  at  z  —  0; 
oz 

.  .  d6\  ,  ,  dd2 

A'l(r)a7  =  A'2<r)a7  at  z  =  “; 

61=62  at  z  =  a; 

62  =  0  at  infinity. 


The  first  boundary  condition  (8)  imposes  the  insulating  boundary  condition  at  the  air/material 
interface  which  assumes  that  no  energy  escapes  into  the  air  at  the  air/material  interface. 
This  is  a  reasonable  approximation  for  most  materials  under  consideration,  since  heat  flow 
by  conduction  through  the  substrate  far  exceeds  the  loss  by  radiation  or  convection  at  the 
air/material  interface.  The  second  boundary  condition  (9)  requires  energy  conservation 
for  heat  flow  across  the  interface  between  regions  1  and  2.  Here  we  have  assumed  that 
K\(T) / K2(T)  =  o  is  a  constant  independent  of  temperature.  This  approximation  can  be 
met  by  many  material  combinations  (e.g.  silicon  on  sapphire).  The  third  boundary  con¬ 
dition  (10)  corresponds  to  temperature  continuity  at  the  interface  between  two  regions, 
and  the  last  boundary  condition  (11)  implies  that  the  temperature  rise  becomes  zero  at 
infinity,  or  equivalently,  the  temperature  of  the  solid  approaches  the  ambient  temperature 
at  infinity. 

The  energy  absorbed  in  the  top  surface  of  region  1  (z  =  0)  is 


(12) 


Q(x,  y,  z,  t)  = 


3Qo 


7rrn 


exp 


/  2? ft  \o  /  7  . 

lx  —  a\  cos - )  +  (y  +  b  1  sin - ) 

_ s * _ s* 

ro/3 


5(z), 


in  the  case  of  a  rotating  laser  beam.  Here  Q0  is  the  power  absorbed  and  it  is  the  product 
of  the  absorptivity  and  the  laser  power  [13] ,  r0  is  the  characteristic  beam  radius  typically 
defined  as  the  radius  at  which  the  intensity  of  the  laser  beam  drops  to  5%  of  the  maximum 
intensity,  s*  is  the  rotating  period.  The  Dirac  delta  function  in  2  expresses  the  assumption 
that  all  the  energy  is  absorbed  at  the  surface.  (12)  describes  a  clockwise  rotating  beam 
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in  the  x-y  plane.  If  cp  =  bi,  then  a\  denotes  the  beam  radius  of  the  rotation  along  the 
z-axis.  When  ai  =  0,  (12)  depicts  the  case  of  a  dithering  laser  beam. 

The  general  solution  of  (6)  and  (7)  can  be  obtained  by  a  Green’s  function  method  [2,  4], 
Denote  the  Green’s  functions  for  regions  1  and  2  as  G\  and  G2 ,  respectively.  Then  the 
Green’s  functions  G'i  and  G2  satisfy 

f)G  -> 

(13)  -  DiAGi  =  —8(r  —  r')5(t  -  t), 
and 

(14)  -  D2AG2  =  0, 

where  | f  —  r'\  —  (x  —  x)2  +  (y  —  y')2  +  (z  —  z')2.  The  initial  conditions  are 

(15)  G'i  =  0  at  t  —  t  for  r^r';  6^2  =  0  at  t  =  t, 


whereas  the  boundary  conditions  for  G'i  and  G2  are  the  same  as  (8)- (9)  with  6\  replaced 
by  G'i  and  02  replaced  by  G2. 

Following  the  work  of  [2,  4],  one  can  express  the  Laplace  transformed  Green’s  functions 


as 


(16) 


G\  — 


G2  = 


2nDi 


2nDi 


exp  (—rjiz)  + 

exp(-r72^)  + 


(arii  ~  rte)  exp(— r)\d)  cosh(riiz) 


(Mm) 


arj i  sinh(Tyia)  +  r\2  cosh^a)  J  iji 

(arii  ~  V2)  exp(-ri2z)  cosh(?]i a)" 
arii  sinh(pia)  +  r/2  cosh(?]ia) 


r  (  m^o(^-R)  ac 

x  exp  [a  (772  -  7]i)J — - - di 


m 

where  G 1  and  G2  are  the  Laplace  transforms  of  G\  and  G2 ,  respectively: 
(17) 

and  Jo  is  the  Bessel  function  of  tbe  first  kind  or  order  zero, 

Vi  =  \[i2  +  n7,  *  =  1)  2, 


Gi=  exp[— p(t  —  t')]Gid(t  —  t'),  i  —  1, 2, 

Jo 


d£ 


(18) 


R  =  y/{x-  xdy  7  (//  i/y, 


a  =  K\(T)  /  K2(T)  =  constant. 
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In  the  appendix  we  verify  that  the  expressions  of  G\  and  G2  in  (16)  are  correct. 

With  the  help  of  the  three-dimensional  Green’s  function,  the  solution  of  the  inhomoge¬ 
neous  heat  equations  (6)  and  (7)  for  a  semi-infinite  medium  can  be  expressed  as 


91  (x,y,z,t)  = 


27r2r02H1iG(To)  J 


rt  ,0O  ,0O  ,0O  rS+iOQ  exp^(t  _  g')] 

— oo  J— oo  J — oo  Jo  J 5— zoo  27fi 


x  "exp(— ?7i^)  +  ~  ^2)  exp(-?71a)  coshfaz)'  £J0{£y/(x  -  x')2  +  (y  -  y')2) 

(19)  L  1  sinh(r^ia)  +  ?72  cosh(r^ia)  J  77 1 

/  ^  27 vt  .0  /  <  ,  .  2nt  .0 

(a;  —  ax  cos  — )  +  (y  +h  sin— J  ;  ;  , 

x  exp - - - 5- - - -  dp  d£  dy  dx  dt 

ro/3 


for  region  1  and 


(20) 

92(x,y,z,t )  = 


3Qo 

27r2r2D1^2(T0) 


rt  r  r  r  rs+ioc  exP[p(t  -  Q] 
—00  J — 00  a — 00  ■ '  ft  J  A — 27T7 


(«77i  -  ?72)exp(-772^)cosh(?71a)l  ^1^o(Cv/(^--^7F+Td--¥F) 

x  exp(— 772^)  + - — j- - 7- - rr - y-  exp  a  772 -771)  - - - 

«77i  smh(77ia)  +  772  cosh  (771  a)  J  771 


/  1  2irt  >2  /  '  ,  .  27t t  .  2 

(x  —  ai  cos  — — )  +  (y  +  61  sm  — -) 


x  exp 


dp  dt;  dy  dx  dt 


for  region  2. 

The  formulas  in  both  (19)  and  (20)  involve  quintuple  integrals.  In  order  to  reduce 
them  to  expressions  with  double  integrals,  we  apply  the  same  procedure  as  in  [4],  First, 
we  consider  the  integral  over  x  and  y  in  (19)  which  is  denoted  as  Ip. 
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Introducing  the  new  independent  variables  X  =  x  —  x  and  Y  —  y  —  y  into  (21),  we  obtain 


Jo  (VUUTW) 


(X  —  x  +  ai  cos - )  +  (y  —  y  —  b\  sin - ) 


x  exp 


dY  dX 


Jo  (VVY2  +  F2)  exp 


I -Rn.fi  —  LI 


JF  dX. 


—>  2nt  —>  2nt  —>  —>  —>—> 

Here  Rab  =  (x  —  0  \  cos  — —  )i  +  (y  +  ly  sin  — —)j  and  L  =  Xi  +  Y j .  It  is  convenient  to 


write  (22)  in  polar  coordinates: 


J0  (^L)  exp  - 


Rnh  +  L2  -  2 RabL  cos  6 


Ldcf)  dL 


2 RnhL , 


=  2n  Jo(£L)Iq(  J^r)exp  — 


Rlh  +  L2 


L  dL, 


2j[t  2i vt 

where  R2fe  =  (x  —  oq  cos - )2  +  (y  +  b\  sin - )2,  L2  =  X2  +  Y2  and  Jo  is  the  modified 

s  s 

Bessel  function  of  the  first  kind  defined  by  Io(z)  =  4  C  exp[z  cos  6]d0. 

In  order  to  reduce  (23)  further,  we  employ  an  identity  from  [10]  (Page  963,  Eq.  (8.5.10) 
and  Eq.  (a)): 


Lexp 


-  .+2X°  Lj(tpw)  =  f  J0(kx)exp[-p2k2]J0(kx0)kdk. 

4 P2  J  2pz  Jo 

-  f?  ,  onrl  rrr^  -  nr*2  /I  o  i  to  frv  (O  /U  TTrn  oKf om 


Introducing  x  =  L,  x0  =  i?a6  and  p2  =  Tq/12  into  (24),  we  obtain 

T  2  _i_  d2  1  O  TD  T  r2  poo 

(25)  exp - Iq(-*)  =  JL  J0(kL)  exp[— r^k2  / 12]  J0(kRab)kdk. 

r0/6  J  r0/u  o  Jo 


Inserting  (25)  back  in  (23)  ,we  find 


J  =nrl  r  r 


Jo(£L)J0(kL)J0(kRab)k  exp 


Pr2l 

—2  LdLdk. 


Substituting  the  orthogonality  property  of  Bessel  functions  (also  called  the  closure  equa¬ 
tion) 

(27)  j™  LJ0(kL)J0(£L)dL  =  -2=d(k  -0  =  l$(k  -  £), 
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in  (26)  yields 


(28) 


/,  = 


7rrn 


5(k  -  £)J0(kRab)exp 


k2rl 

'TP 


71V 

dk  =  Jo{^Rab)  exp 


12 


We  may  now  return  to  (19).  The  linear  heat  rise  in  region  1  can  be  written  as 

Q0  V  f°°  f5+i°° 


91  (x,y,z,t)  = 


(29) 


x 


cxp(  rftz)  + 


4ir2iDiKi(T0)  J_00  J0  J5_ioo 

(ary  —  rfe)  exp(— ?7ia)  coshOiz) 


exP [p(t  -  t )] 


£  Jo(^R-ab) 
x - exp 

Vi 


ary  sinh(r7ia)  +  772  cosh(?7ia) 
dp  d£  dt 


erf 

12 


In  [4]  a  scanning  laser  beam  was  used  and  the  steady  state  solution  was  studied.  With  the 
help  of  the  residue  theorem  one  was  able  to  find  the  integral  with  respect  to  the  variable 
p  and  thereby  reduce  the  triple  integral  to  a  double  integral.  Here  a  rotating  or  dithering 
laser  beam  is  considered  and  we  are  interested  in  time-dependent  solution.  As  a  result, 
the  techniques  used  in  [4]  cannot  be  extended  here. 

Now  making  a  change  of  variables  t  —  t  —  t' ,  we  can  rewrite  (29)  as 

Qo 


poo  poo  p5-\-ioo 

/  /  /  exp  [pt\ 

J  0  Jo  J  5—ioo 


9i(x,y,z,t)  = 


X 


(30) 


exp(-r/iz)  + 


An2iDiKi(T0) 

(aip  —  772)  exp(— ?7ia)  cosh^iz) 


^Jo(^Rab) 
x - exp 

Vi 

Qo 


ary  sinh(?7ia)  +  r/2  cosh  (771  a)  J 
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where  =  (x  —  ai  cos  — —  )2  +  (y  +  bi  sin - - - )2,  Fi(i,  £)  is  the  inverse  Laplace 


s *  s * 

transform  of  a  complicated  function  given  explicitly  by 

p(5+zoo 


(31)  Fi(i,0  = 


1 

2tTX  J  £ 5—ioo 


ex 


P  \pt\ 


exp  ( -7/1  z)  + 


(«77i  —  772)  exp(— 771a)  cosh^iz) 


api  sinh(77ia)  +  772  cosh  (771  a) 

and  recall  that  ry  =  p£2  +  p/ Dt,  i  =  1,2. 

In  a  similar  fashion,  the  temperature  for  region  2  can  be  expressed  in  the  form 


m 


-dp, 
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92(x,y,z,t)  = 
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POO  POO 


2nDiK2(T0) 
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'0  Jo 


pp 

12 
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where  F2(i,£)  is  the  inverse  Laplace  transform 


(33) 


F2(i,0 


exp(-r]2z)  + 


(/yr/i  ~  V2)  exp(-ri2z)  cosh(r]ia) 
ar/i  sinh(r]ia)  +  r/2  cosh^a) 


x  expjafe  -m)] 
Vi 


The  analytical  solutions  above  involve  the  inverse  Laplace  transform  and  multiple  integrals 
and  thus  it  is  numerically  difficult  to  obtain  meaningful  solutions.  As  a  future  work  we 
will  try  to  simplify  the  expressions.  In  the  next  section  we  will  solve  the  heat  equations 
directly  in  a  finite  domain. 

3.  Numerical  solutions  in  a  finite  domain 


In  this  section  we  use  the  commerical  software  COMSOL  to  solve  the  heat  equations 
directly  for  a  two-layer  structure  in  a  finite  domain.  In  the  simulations  below  we  model  a 
two-layer  structure  where  the  substrate  is  made  of  aluminum  and  the  thin  coating  layer 
is  made  of  copper. 

Figure  2  depicts  the  model  geometry  for  a  two-layer  structure.  It  consists  of  a  substrate 
made  of  one  material  and  a  thin  layer  of  coating  of  different  material.  Ideally  the  coating 
materials  should  have  larger  thermal  conductivity  than  the  substrate  materials.  We  as¬ 
sume  that  there  is  no  energy  escaping  into  the  ambient  at  the  air/structure  interface  so 
the  boundaries  are  treated  as  insulating  boundaries.  For  most  situations  this  is  a  reason¬ 
able  assumption  since  heat  loss  by  radiation  or  convection  at  the  interface  is  often  much 
smaller  than  heat  flow  by  conduction  through  the  material.  Furthermore,  the  structure  is 
assumed  to  have  the  same  temperature  as  the  ambient  initially  so  the  initial  temperature 
rise  is  set  to  zero. 

For  all  the  computations  presented  in  this  paper,  we  choose  aluminum  as  the  substrate 
material  and  copper  as  the  coating  material.  Their  thermal  properties  are  listed  in  Table 
1.  For  laser  beams  described  in  (12)  we  choose  Qo  =  1.0  x  10 5W/m2,  r0/V6  =  0.02 m, 
s*  =  Is,  cq  =  0.25m,  b\  =  0m  (for  a  dithering  laser  beam),  b\  =  0.25m  (for  a  rotating  laser 
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laser  beam 

V 


Figure  2.  A  plot  of  a  two-dimensional  two-layer  structure. 

beam).  The  two-layer  structure  is  a  three-dimensional  box  with  dimension  1  m  x  lm  x  lm 
and  the  center  of  the  laser  beam  hits  the  center  of  the  top  surface. 


Table  1:  Thermal  properties  of  copper  and  aluminum 


Property  Name 

Copper 

Aluminum 

specific  heat  (J / (kg  ■  K)) 

385 

900 

thermal  conductivity  (W/(m-  K)) 

400 

160 

thermal  diffusivity  (m2/s) 

1.1942  x  10-4 

6.5844  x  10“5 

melting  point  ( K ) 

1356 

933 

density  (kg/m3) 

8700 

2700 

In  Figure  3(a)  we  plot  the  temperature  rise  in  an  aluminum  box  induced  by  a  Gaussian 
beam  at  time  =  Is  whereas  Figure  3(b)  shows  the  temperature  rise  on  a  horizontal  slice 
of  Figure  3(a)  through  z  =  0.99 m.  The  maximum  temperature  ( 1609 AT)  is  reached  at  the 
top  surface  where  the  beam  hits  the  box.  The  maximum  temperature  rise  drops  to  937 K 
at  the  horizontal  slice  0.01m  below  the  top  surface.  Figure  3(c)  shows  the  temperature 
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rise  in  an  aluminum  box  coated  with  0.1%  thin  copper  layer  on  the  top  surface.  The 
maximum  temperature  rise  falls  to  13527%  which  is  about  16%  reduction  from  the  case 
without  coating.  Figure  3(d)  depictes  the  corresponding  temperature  rise  at  the  horizontal 
slice  0.01m  below  the  top  surface.  Now  the  maximum  temperature  is  reduced  to  797 K. 


Figure  3.  (a)  Temperature  rise  of  an  aluminum  box  induced  by  a  Gaussian 
beam  with  no  coating;  (b)  a  horizontal  slice  of  part  (a)  through  z  =  0.99; 
(c)  Temperature  rise  of  an  aluminum  box  with  0.1%  copper  skin  on  the  top 
surface;  (d)  a  horizontal  slice  of  part  (c)  through  z  =  0.99 
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Figure  4(a)  shows  further  the  dependence  of  the  maximum  temperature  rise  of  the  two- 
layer  box  on  the  thickness  of  the  coating  layer.  As  expected,  the  maximum  temperature 
rise  decreases  when  the  thickness  of  the  coating  layer  increases.  Figure  4(b)  predicts 
the  asymptotic  behavior  of  the  maximum  temperature  rise  (in  K)  as  a  function  of  the 
thickness  of  the  coating  layer  (in  mm)  using  the  least  squares  fitting: 


(34) 


T  = 

max 


-3290.32  3459.76  _ 

- a - + - +  252.79. 

az  a 


(a)  (b) 

Figure  4.  (a) The  maximum  temperature  rise  of  the  two-layer  aluminum 
box  induced  by  a  Gaussian  beam  as  a  function  of  the  thickness  of  the  copper 
coating  layer.  The  computed  data  are  represented  by  the  red  diamonds 
whereas  the  solid  curve  is  obtained  by  a  cubic  spline  interpolation.  (b)The 
asymptoic  behavior  of  the  maximum  temperature  temperature  rise  ( Tmax ) 
as  a  function  of  the  thickness  of  the  coating  layer  (a).  The  solid  curve 
corresponds  to  the  least  squares  fitting  function  Tmax  =  ~32920,32  +  3459,76  + 

252.79. 

Next  we  let  the  laser  beam  dither  or  rotate  and  study  the  effect  on  the  maximum 
temperature  rise.  We  first  apply  a  dithering  Gaussian  beam  on  the  top  surface  of  the 
aluminum  box  without  coating  and  the  temperature  rise  is  depicted  in  Figure  5(a).  Then 
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we  apply  the  same  dithering  Gaussian  beam  on  the  top  surface  of  the  aluminum  box 
coated  with  0.1%  thin  copper  layer  on  the  top  surface.  The  corresponding  temperature 
rise  is  shown  in  Figure  5(b).  Due  to  the  effect  of  the  dithering  beam  the  maximum 
temperature  rise  falls  to  329 K  and  it  is  further  reduced  to  268 K  if  a  thin  layer  (0.1%) 
of  coating  is  applied.  For  the  model  case  considered  here,  moving  beam  is  more  efficient 
than  coating  in  reducing  maximum  temperature  rise. 


Time=l  Surface:  Temperature  (K) 


A  329.4 


Time=l  Surface:  Temperature  (K) 


250 


150 


100 


▼  -3.111 


(a)  (b) 

Figure  5.  (a)  The  maximum  temperature  rise  of  an  aluminum  box  induced 
by  a  dithering  Gaussian  beam.  (b)The  maximum  temperature  rise  due  to 
a  dithering  Gaussian  beam  on  an  aluminum  box  coated  with  0.1%  copper 
layer. 


Figure  6  is  similar  to  Figure  5,  but  for  a  rotating  Gaussian  beam.  The  maximum 
temperature  rise  induced  by  a  rotating  beam  is  smaller  than  the  one  induced  by  a  dithering 
beam. 

The  dependence  of  the  maximum  temperature  rise  on  the  thickness  of  the  coating 
layer  induced  by  a  dithering  or  rotating  laser  beam  is  demonstrated  in  Figure  7(a)  and 
(b)  respectively  together  with  least  squares  fitting  curves.  Figure  7(a)  shows  that  for  a 
dithering  beam  the  maximum  temperature  rise  depends  on  the  thickness  of  the  coating 
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(a) 


(b) 


Figure  6.  (a)  The  maximum  temperature  rise  of  an  aluminum  box  induced 
by  a  rotating  Gaussian  beam.  (b)The  maximum  temperature  rise  due  to  a 
rotating  Gaussian  beam  on  an  aluminum  box  coated  with  0.1%  copper  layer. 


layer  as  a  function 


(35) 


T  = 

^  max 


-15.5595  72.7360 

- - - t - b  210.9354, 


cr 


a 


whereas  for  a  rotating  beam  the  dependence  is  different: 


(36) 


T  = 

max 


-5.9499 


az 


28.0879 

a 


64.0152, 


as  shown  in  Figure  7(b). 

Finally,  for  comparison  purpose,  we  summarize  the  results  of  Figures  3,  5  and  6  in 
Table  2  which  shows  quantitatively  how  the  maximum  temperature  rise  can  be  reduced 
for  various  cases.  From  Table  2,  we  summarize  our  observations:  (1)  rotating  an  object 
reduces  the  maximum  temperature  rise  more  than  dithering  the  object;  (2)  moving  an 
object  is  more  efficient  than  coating  in  reducing  maximum  temperature  rise.  Finally,  it 
is  obvious  from  Table  2  that  the  most  effective  way  to  reduce  the  maximum  temperature 
rise  is  to  combine  coating  with  a  rotating  Gaussian  beam.  From  a  practical  defense  point 
of  view,  this  implies  that  in  order  to  reduce  the  maximum  temperature  rise  in  an  object 
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(a) 


(b) 


Figure  7.  The  maximum  temperature  rise  of  the  two- layer  aluminum  box 
as  a  function  of  the  thickness  of  the  copper  coating  layer.  The  computed 
data  are  represented  by  symbols  whereas  the  solid  curves  are  obtained  by 
the  least  squares  approach,  (a)  The  beam  is  a  dithering  Gaussian  beam  and 
the  fitting  function  is  Tmax  =  ~15a2595  +  <2'^360  +  210.9354.  (b)  The  beam  is 
a  rotating  Gaussian  beam  and  the  fitting  function  is  Tmax  =  ~5a9499  + 28,98,9  + 
64.0152. 


hit  by  a  laser  beam,  one  possible  way  to  protect  the  object  is  to  apply  a  thin  layer  of 
coating  materials  on  the  object  and  also  rotate  the  object. 


Table  2:  Summary  of  maximum  temperature  rise  in  various  cases 


Gaussian  beam 

Dithering  Gaussian  beam 

Rotating  Gaussian  beam 

rri no  coating 
-L  max 

1609 

329 

109 

( K ) 

(reference  temperature) 

(79.55%  reduction) 

(93.23%  reduction) 

rp 0.1%  coating 
^max 

1352 

268 

86 

( K ) 

(15.97%  reduction) 

(83.34%  reduction) 

(94.66%  reduction) 
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4.  Conclusions 

In  this  paper  we  have  derived  the  temperature  rise  induced  by  a  rotating  or  dithering 
Gaussian  laser  beam  for  a  semi-infinite  two-layer  structure.  Numerical  solutions  for  a 
finite  two-layer  structure  are  provided  using  COMSOL.  It  is  found  that  the  maximum 
temperature  rise  can  be  reduced  significantly  by  applying  a  thin  layer  of  coating  or  by 
moving  the  beam  (or  object). 
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Appendix.  Verification  of  the  Green  Functions 


In  this  appendix  we  verify  that  the  Lapace  tranformed  Green  functions  in  (16)  are 
correct.  For  simplicity  we  consider  only  the  case  r  =  0  (otherwise  one  can  make  a  shift 
to  obtain  this  case). 

First,  we  check  the  expression  of  G\.  Taking  the  Laplace  transform  of  (13)  and  using 
the  initial  condition  (15),  we  have  an  equation  for  Gp 

(37)  pGi  —  DiAG\  =  8(r  —  r). 


Using  the  cylindrical  coordiates,  we  can  express  (37)  explicitly  as 

d2Gx 


(38) 


pG  i  —  D\ 


d2G1  1  dGl 

+ 


=  m 


dR 2  RdR  dz 2 

We  consider  the  case  where  r  ^  0.  For  G i  given  in  (16),  it  is  straightforward  to  obtain 
that 


(39) 
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^  n  ( d2Gi  ,  19G,  ,  a2G, 
W+W+  ~SF 


f  f  \  ,  ( arh  -  m)  exp(-77ia)  cosh  (771  z)  £ 

-  /  exp(—r}iz)  H - —7 - 7 - T7 - —  — 

r  Jo  L  arji  sinh[r)ia)  +  r)2  cosh{r)ia)  J  ip 

d2Jom  ,  ,, 


R  dR 


+euiR)  di 


Introducing  a  new  variable  s  =  £ R ,  we  get 


j0  m  =  Jo(s), 


dJ0  _  dJ0  OR  _  1  dJ0 
~ds  ~  ~dR~ds  ~  l~dR 

d2j0  _  1  a2  J0 
Gh2  ~  ^aR2" 


It  follows  immediately  that 

,,,,  guo m  1 9Jom 

[  ’  am  r  dR 


+  £2  Jo{£R)  —  Jp  s2^0  (s)  +  sJq(s)  +  s2Jo(s)  —  0 


where  in  the  last  step  we  have  used  the  fact  that  J0  is  the  Bessel  equation  of  order  zero. 
Substituing  (42)  in  Eq.  (40),  we  obtain,  finally, 

(43)  pGi  —  D1AG1  =  0  for  r  7^  0. 

In  a  similar  way  one  can  confirm  that  G2  satisfies 

(44)  pG2  -  D2AG2  =  0. 

Finally,  we  check  the  boundary  conditions. 

OGi  1  f°°  [  .  .  (0:771  -  p2)  exp(-r^ia)  sinh^z)]  A  T  IA 

(45)  ~R~  =  R—FT  \  -exp(-?7i  z)  +  - - — - - — - — - €J0(€R)d£ 

oz  2nL>i  J0  |_  arh  smh(?7ia)  +  rj2  cosh(?7ia) 


This  gives 


2  poo  2  r°° 

sMs)ds, 


where  we  have  made  a  change  of  variables  s  =  ^R.  Since  Jq(s)  satisfies  Bessel’s  equation 


of  order  zero,  we  have 


(47) 


S2Jq  (5)  +  sJq  +  52Jq(5)  —  0 
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or, 

(48) 

Applying  this  result  to 

(49) 


SJ„(S)  =  -K(S)  +  j;(S)]  =  -^K(S)]. 


ggl 

dz 


z=0  from  Eq.  (46),  we  find  that 
0Gl 


dz 


z=0  —  0, 


provided  that  sj'0(s )  — >  0  as  s  — »  oo. 
Note  that  from  (16),  we  obtain 
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Consequently,  (?i|2=a  =  G2\z=a. 

Now  we  take  the  partial  derivative  of  G2  with  respect  to  2,  thereby  obtaining 
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On  the  other  hand,  from  (45)  we  obtain 

(ar/ 1  —  172)  sinh(?7ia) 
ar/i  sinh(?7ia)  +  r]2  cosh(?7ia) 

It  then  follows  immediately  that 
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or, 

fril  K  (rp\dG^\  _  K  frp\dG2 

(54)  Ki (T)  ^  \z=a  —  K2(T)  I z=a- 

recalling  that  a  =  Ki(T)  /  K2(T) .  The  boundary  condition  Go  -A  0  at  inhnity  (i.e. 
z  — »  00)  is  obvious. 


